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The role of sound in the dynamics of mesoscale systems is typically ne- 
glected, since frequently the associated time scales are much smaller than 
all the other time scales of interest. However, for sufficiently small objects 
embedded in a solvent with a sufficiently small sound velocity, sound can 
play a crucial role. In particular, behavior resembling viscoelasticity has 
been theoretically predicted for non-viscoelastic fluids. This effect is due 
to the interference of the propagation of sound waves caused by the solute 
particle's motion and hydrodynamic vortex formation. We demonstrate 
this effect, known as backtracking, in computer simulations employing the 
method of multiparticle collision dynamics. We systematically study the 
influence of sound on the dynamics of the solute particle, and find that it 
disappears in the long-time limit. Thus, we confirm that sonic effects at 
the single-particle level can be neglected at sufficiently long times. 

1 Introduction 

The progress in single molecule tracking techniques provides a route to study 
the dynamical aspects of individual particles, such as colloids, polymers, vesicles 
or even cells. Specifically, insight is achieved in the influence of the environment 
on the particle behavior. This is of particular importance for dilute suspen- 
sions, where hydrodynamic effects determine the transport properties of the 
solute particles. Here, the solvent characteristics, e.g., compressibility, are of 
paramount importance. 

Investigations of the dynamics of small objects suspended in a medium and 
driven by thermal fluctuations, i.e. undergoing Brownian motion, have a long 
fruitful history [T|. In the classical picture of Brownian motion, the dynamics 
of a solute particle experiencing viscous drag and random collisions with the 
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Setup 


Box size 


Collision angle 


Solvent density 


Collision time 


Colloid radius 


Colloid mass 




L [a] 


a [deg.] 


Pf [raja^] 


5t [y/ma'^/kBT] 


R [a] 


M [to] 


1 


32 


130 


10 


0.01, 0.02, 0.05, 0.1 


2.0 


335.1 


2 


32 


90 


5 


0.005, 0.01, 0.05 


2.0 


125 



Table 1: Parameters of the different setups for which simulations have been 
performed. This choice of parameters allows to illustrate the smooth transition 
from normal behavior to backtracking. 



molecules of the solvent is described by the Langevin equation, 

Mi(t) = -7i(<) + r(<). (1) 

Here M is the particle mass, x its coordinate, 7 is the solvent's viscous friction 
coefficient and T{t) is a stochastic force which models collisions of the solute 
particle with the molecules of the solvent. The stochastic force is assumed to 
be stationary, Markovian, and Gaussian with zero mean (white noise), i.e. 

{m) - 0, (2) 

{T{t)V{t')) = 2kBT-f5{t - t') . (3) 



Setup 


Collision time 


Viscosity r] 


Time scale 




St \^m.a^lkBT\ 


[kBTm/a'^] 


ratio Tf/Ts 


1 


0.01 


82.2 


0.3 




0.02 


41.2 


0.6 




0.05 


16.7 


1.5 




0.1 


8.7 


3.0 


2 


0.005 


44.5 


0.3 




0.01 


22.3 


0.6 




0.05 


4.6 


2.8 



Table 2: Viscosities and ratio of the viscous to the sonic time scales correspond- 
ing to the different collision times of the two setups. 

The Langevin equation ([T]) is a stochastic version of Newton's equations 
of motion with an effective force which models random uncorrelated collisions 
between the solute particle and the solvent molecules. 

It is well known that this simple picture of Brownian motion is largely in- 
complete for objects in dilute solution, in particular for short times. As early as 
in the 1960's, Alder and Wainwright [2 [3] observed by means of computer sim- 
ulations strong deviations of the behavior of a Brownian particle from behavior 
predicted by the simple Langevin equation ([T]). These deviations arise from 
hydrodynamic correlations. Furthermore, Zwanzig and Bixon |3] showed that 
in order to satisfy the equipartition theorem, it is necessary to account for the 



2 



finite compressibility of the medium. Since then, correlated momentum trans- 
port both by shear and by longitudinal wave propagation is known to play an 
important role in fluid dynamics. [5]. Assuming fully-developed hydrodynamics, 
Stokes derived for a sphere in an incompressible fluid 7 = GnriR for stick bound- 
ary conditions at the solute-solvent interface, or 7 = An-qR for slip boundary 
conditions in the long-time limit, based on a time-dependent response function, 
which relates the friction coefficient with the size of the colloid [B]. Here, R is 
the radius of the sphere and is the shear viscosity of the solvent. 

However, the Stokes friction coefficient describes the friction experienced by 
the colloidal particle only at long observation times. At short times, viscous 
and compressibility effects strongly manifest themselves. This implies that at 
different stages of its motion, the solute particle experiences friction of different 
magnitude. It turns out possible to calculate the resulting friction coefficient in 
Fourier space: the corresponding frequency-dependent expressions taking into 
account viscous vortex formation were first provided by Stokes [51 [7] , and since 
then refined to include compressibility, an arbitrary degree of slip [H [HI (THl E] 
and several geometries [T^ . 

While viscous effects have received broad attention in the last years, sonic 
effects in viscous fluids are usually neglected. This can be easily understood in 
terms of the relevant time scales [IS]: 

• the sonic time scale corresponding to the time it takes a sound wave to 
propagate one colloid radius, = R/v^, with Vg the velocity of sound in 
the solvent; for a micrometre-sized colloid in water 1 ns; 

• the inertial time scale corresponding to the characteristic time of ballistic 
motion of the colloid, tb = Af/7; for the same colloid in water tb ~ 0.1 

• the viscous time scale corresponding to the time of momentum diffusion 
by one colloidal radius, t/ = R? jv^ with v the kinematic viscosity of the 
solvent, V = rj/ pf, where /O/ is the mass density of the solvent; in the same 
system as above t/ « 1 /is. 

Thus, in conditions typical in a soft-matter experiment the sonic time scale is 
many orders of magnitude smaller than the other time scales of interest, and it 
is therefore assumed that sonic effects can be neglected. 

However, recent experimental investigations have highlighted the role of 
sound in the formation of hydrodynamic correlations in many-body systems |14( 
[TUl [T5] . Previously, it was assumed that cross-correlations between pairs of col- 
loidal particles would build up on the viscous time scale r/, and only recently 
it has been shown that instead they build up on the sonic time scale Tj. |15] . 
Computer simulations revealed that the origin of this phenomenon lies in mul- 
tiple scattering of sound waves [16] , implying that the role of sound effectively 
extends to times much longer than the sonic time. Such phenomena highlight 
the importance of a solid theoretical understanding of the interplay between the 
sonic and the viscous effects. 
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Figure 1: Schematic representation of MFC. (A) Each solvent particle undergoes 
streaming, whereby its coordinates are updated according to its velocity and 
the collision time. (B) Farticles are sorted into collision cells. Within each cell, 
the velocity of each particle is rotated with respect to the cell's centre-of-mass 
velocity by a fixed angle around a randomly oriented axis. (C) The structure of 
the collision grid is shifted to guarantee Galilean invariance. (D) All particles 
which may interact with the solute object do not undergo streaming. Instead, 
their equations of motion are integrated with a much finer molecular dynamics 
time step. 



Even more recent theoretical investigations have revealed other ways in 
which sonic and viscous effects can interfere [17] . demonstrating that sound 
propagation can give rise to non-trivial effects not only at a finite concentration 
of solute particles, but also at the single-particle level. More specifically, it was 
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Figure 2: Left: computer simulations for Setup 1. The solid line is a guide 
for the eye showing the t^'^l'^ slope of the long-time tail of the velocity auto- 
correlation function which indicates that hydrodynamic correlations have fully 
developed. The dashed line in the inset shows the Enskog exponential decay. 
Right: theoretical calculation for Setup 1. Solid lines correspond to calculations 
including viscous effects and solvent compressibility. Dotted lines show the cal- 
culations in the incompressible limit. The dashed line in the inset shows the 
Enskog exponential decay. 



shown that under conditions where the compressibility of the fluid is sufficiently 
large and its bulk viscosity sufficiently small, the generation of sound waves 
during the motion of a single colloidal particle may lead to a hydrodynamic re- 
action force which can strongly influence its future motion. This effect, known 
as backtracking |17j , is manifested in the velocity autocorrelation function of the 
colloid. Instead of a monotonic decay, it exhibits a kink at relatively small sol- 
vent compressibility, or even crosses zero at high solvent compressibility. Thus, 
it is predicted that the colloid should change its initial direction of movement at 
times corresponding to the viscous and the sonic times, which should be of the 
same order of magnitude. The backtracking effect is reminiscent of viscoelastic- 
ity, but it has a very different, dynamical origin, which can be mathematically 
described as follows. 

The decay of the velocity autocorrelation function C„t,(<) of a solute particle 
is given by 

a,(t) = fcBr$(t) (4) 

with the response function 

m = ^ J yi^)e-''^'di., (5) 

— oo 

where y{co) is the admittance tensor, 

and the frequency-dependent friction coefficient j{uj) can be calculated using 
the analytical formulae of Refs. [8l[12]. Essentially, it was demonstrated [17] 
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Figure 3: Left: computer simulations for Setup 2. The solid line is a guide 
for the eye showing the t^'^l'^ slope of the long-time tail of the velocity auto- 
correlation function which indicates that hydrodynamic correlations have fully 
developed. The dashed line in the inset shows the Enskog exponential decay. 
Right: theoretical calculation for Setup 2. Solid lines correspond to calculations 
including viscous effects and solvent compressibility. Dotted lines show the cal- 
culations in the incompressible limit. The dashed line in the inset shows the 
Enskog exponential decay. 



that, provided that y(y>) exhibits certain analytical properties in the complex 
plane, velocity reversal of the solute particle will occur. It was shown that this 
velocity reversal is due to the elasticity of the fluid corresponding to the finite 
velocity of sound. 

In this contribution, we address the question to what extent sonic effects con- 
tribute to the dynamics of colloids on the single-particle level. To this end, we 
study backtracking in computer simulations which capture thermal fluctuations, 
hydrodynamic vorticity and sound propagation. We investigate the influence of 
sound on long-time observables, namely the diffusion coefRcient of a colloidal 
particle suspended in a solvent. The results are timely, since time scales in 
the range of the typical hydrodynamic time scales in soft-matter systems - typ- 
ically of the order of microseconds - are becoming more and more accessible 
experimentally. 



2 Simulation Method 

We employ the method of multiparticle collision dynamics (MFC) to investi- 
gate backtracking and its influence on the long-time observables in a simple 
system of a spherical colloid embedded in a viscous solvent. MFC is known 
to correctly capture both thermal fluctuations and hydrodynamics on coarse- 
grained scales [IHl im HO] , has an H-theorem and satisfles the fluctuation 
theorem [22]. MFC is ideally suited for observation of backtracking since it 
satisfies the necessary criteria p^: the bulk viscosity is zero and the sound 
velocity Vg is small, since the MFC fluid behaves as an ideal gas. Therefore, 
the simulations present ideal conditions for the observation of the interference 
of sonic and viscous effects, since the sonic, viscous and inertial time scales are 
very similar. 
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To model the solvent in the simulations, a cubic box of lateral size L with 
periodic boundary conditions is filled with N pointlike solvent particles of mass 
TO with continuous coordinates fi and velocities Vi, i G [1;^]- The average 
kinetic energy per solvent particle determines the thermal energy of the solvent. 
At each step of the simulation, the solvent undergoes streaming and collisions. 
During streaming, the coordinates of each solvent particle are updated according 
to its velocity Vi and the collision time St: — + ViSt ([l] (A)). Then, all 
particles are sorted into collision cells of lateral size a according to their positions 
within the simulation box ([l] (B)). The collision cells form a cubic collision grid 
spanning the whole simulation box; in order to guarantee Galilean invariance, 
the structure of this collision grid is shifted at each simulation step [201 US] 
([I] (C)). Within a collision cell, the velocities of the particles are rotated with 
respect to the cell's centre-of-mass velocity by a fixed rotation angle a around a 
randomly oriented axis. Energy and momentum are thus conserved locally. Due 
to the existence of an H-theorem [3T], the velocity distribution of the solvent 
particles relaxes to a Maxwell-Boltzmann distribution in just a few simulation 
steps, regardless of the initial distribution of velocities. 

We couple the solute colloidal particle to the solvent employing the stan- 
dard hybrid molecular dynamics approach |24j : solvent particles interact with 
the solute through a shifted Lennard- Jones potential with a cutoff radius corre- 
sponding to the colloidal radius R. This form of coupling implies slip boundary 
conditions at the solute-solvent interface. Stick boundary conditions can also 
be studied with MPC [551 

We choose the solute-solvent interaction potential to be 



All solvent particles which may interact with the colloid within a MPC step 
do not undergo streaming. Instead, their equations of motion are integrated 
with the smaller molecular dynamics (MD) time step, taking into account the 
solvent-solute interaction potential ([l] (D)). Afterwards, they are sorted into 
collision cells and undergo collision with the rest of the MPC solvent particles. 

We choose the following normalization: solvent particle mass to = 1, thermal 
energy ksT = 1 and collision cell size a = I. Thus, mass is measured in units 
of TO, energy in units of fc^T, length in units of a and time in units of to = 
y/ma^ /kBT. We perform simulations for several setups. The corresponding 
setup-specific parameters are listed in [T] For the observation of backtracking, it 
is important that the sonic time scale is not too small, i.e. the ratio of viscous to 
sonic time scales is not large. In [2] we list the shear viscosities corresponding to 
each of the time steps of the two setups; the viscosities are calculated using the 
analytical expressions of Refs. [IHl dOl HZl UHl UHl [30] • The corresponding ratios 
of the viscous time scale r/ = B?pf /r] to the sonic time scale Ts = R/vg are also 
shown. Since our fluid obeys ideal-gas equations of state, the velocity of sound - 
and, consequently, the sonic time - is fixed, i.e., Vs = ^5/3 [201 [3l]. Therefore, 





7 



xlO' 



„ 10 



Q 5 







• 


Tf/Ts = 0.3 


■ 


Tf/Ts = 0.6 




Tf/Ts = 1.5 




Tf/Ts = 3.0 




1 TTT] — I — rTT| — I — rTT| — I — rTT| — I 

10-2 10-1 10° 10^ 102 
t [ Vma^/keT ] 




— I 1 m~| 1 TTT] 1 TTT"! 1 rTT| 

10-1 10° 10^ 10^ 103 
t [ Vma^/keT ] 



Figure 4: The product of the solvent viscosity rj and the time-dependent diffu- 



sion coefficient D(t) calculated using 10 in simulations (left) and theory (right) 
for Setup 1. Dotted lines show the results of the calculations in the incom- 
pressible limit. In all cases, the long-time limit of the diffusion coefficient is 
unaffected by the effects of sound. 



it is the viscous time scale which changes between the different simulations due 
to the difference in viscosities which, in our case, depend on the collision time 
6t, collision angle a, and the solvent density pf. 

3 Results 

The colloid velocity autocorrelation functions (VACFs), averaged over several 
hundred simulation runs each, are shown in [2] and [3] for Setup 1 and Setup 2, 
respectively. The corresponding theoretical curves are calculated in terms of the 
inverse Fourier transform of the admittance tensor [17] (|5| with the frequency- 
dependent friction coefficient calculated following Ref. [12;. A key ingredient of 
this calculation is the hydrodynamic radius of the colloidal particle R^, such 
that in the zero-frequency limit the Stokes relation is satisfied, ^{Ld = 0) = 
A'ktiRh- Equivalently, the theoretical calculation can be performed if simply 
the long-time limit of the friction coefficient is known. This is the approach 
we adopt: we assume the validity of the Stokes-Einstein relation at long times, 
D = ksT and measure the corresponding friction coefficient 7 in simulations. 

The measured data are in good qualitative agreement with the theory and, 
clearly, several regimes can be distinguished. 

First of all, we observe that at sufficiently long times, all the velocity auto- 
correlation functions exhibit the same long-time tail which depends only on the 
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Figure 5: The product of the solvent viscosity rj and the time-dependent diffu- 



sion coefficient D{t) calculated using 10 in simulations (left) and theory (right) 



for Setup 2. Dotted lines show the results of the calculations in the incom- 
pressible limit. In all cases, the long-time limit of the diffusion coefficient is 
unaffected by the effects of sound. 

properties of the solvent p2l l33 l [34 1. [35 l fT8 ] . 

In this limit, the simulations are in excellent agreement with the theory. We 
again stress that the long-time tail is a direct consequence of the Navier-Stokes 
equation which yields the time-dependent (or, equivalently, the frequency-dependent) 
friction coefficient. The Langevin equation itself predicts a simple exponential 
decay of the VACF. 

At very short times, however, there is an apparent discrepancy between the 
simulations and the theory, as indicated in the insets of|2]and[3j Whereas simu- 
lations exhibit an exponential decay of the velocity autocorrelation function, the 
theory predicts a decay which is faster than exponential. This difference has a 
straightforward explanation. At very short times, the solvent in the simulations 
can not be treated as a continuous medium, and the Navier-Stokes equation can 
not be applied. Instead, one has to assume random collisions between solvent 
particles and the colloid. Here, the Langevin equation applies, and the velocity 
autocorrelation function indeed decays exponentially, with the decay constant 
given by the ratio of the Enskog friction je and the colloid mass M [TS], with 

,.^in-^^feiRK (9) 

3 \ M + m J m 
This behavior is, in principle, a consequence of the coarse-graining approach of 
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Figure 6: Velocity autocorrelation functions for solvent particles for a — 130° 
and p = 10 (left) and a = 90° and p — 5 (right). The points correspond to 
positive values of the VACF and the solid lines correspond to negative values of 
the VACF. 



MFC - the relatively small number of solvent particles and the discrete dynamics 
imply that few collisions occur in a single simulation step, and, consequently, 
a continuum description is not applicable. We remark that the Enskog fric- 
tion coefficient has a different physical origin than the Stokes viscous friction 
coefficient, and that the two are not necessarily related. 

The theory, on the other hand, assumes that a hydrodynamic description is 
valid down to arbitrarily short times. The Enskog regime is completely absent 
from the theoretical treatment. Consequently, the theoretical curves differ only 
when compressibility is accounted for, and only at the sonic time scale (see [2j 
right) . 

We further observe that for characteristic sonic times Tg smaller than viscous 
times Tf, the velocity autocorrelation functions decay monotonically and do not 
exhibit any peculiar properties. However, when becomes larger than Tf, an- 
ticorrelations clearly develop in the velocity autocorrelation function, signaling 
backtracking. The monotonic decay of the VACF is superseded by crossing zero 
and an asymptotic approach of the hydrodynamic t^"^/^ long-time tail. Here, 
the simulation results are in good qualitative agreement with the predictions 
of the theory. We note that if, in the theoretical calculations, the effects of 
solvent compressibility are not taken into account, i.e. if the calculations are 
performed with an infinitely high velocity of sound, no anticorrelations develop 
(dotted lines of [2] and [3|. Therefore, these anticorrelations are indeed a direct 
consequence of the interference of sonic and viscous effects. 

To study the influence of backtracking on long-time transport properties, we 
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calculate the diffusion coefficient as the integral of the velocity autocorrelation 
function, i.e., 

t 

Dit) = J {v{T)v{0))dT . (10) 



Assuming the validity of the Stokes-Einstein relation, 

7 4TTr]R ^ ' 

where D = D{oo), the product Drj depends only on the properties of the colloid, 
namely its size, and does not depend on properties of the solvent, namely its 
shear viscosity, and serves as a good indicator of the relevance of sonic effects 
at long times. 

As shown in |4] and [5] we observe that in the presence and in the absence of 
backtracking, riD{t) tends to the same value in the asymptotic limit. Clearly, 
at the sonic time scale, there are significant deviations between the cases with 
and without backtracking. However, the long-time limit is always the same. 
As before, the simulation results are in good qualitative agreement with theory. 
We note that if the solvent compressibility is eliminated from the theoretical 
calculations, the long-time limit of rjD[t) remains the same (dotted lines of |4] 
and [5]). 

Hence, we have confirmed that even when the effects of sound are very pro- 
nounced and dramatically influence the short-time regime of the solvent parti- 
cle's motion, the net effect remains zero, i.e. the role of sound can be neglected 
when studying long-time observables. 

The main requirement to observe backtracking is that the sonic time scale 
should be equal to or slightly larger than the viscous time scale. This is clearly 
not the case in typical conditions of a micron-sized particle suspended in a liquid 
like water since, as we demonstrated above, in this case the sonic time scale is 
of the order of nanoseconds, whereas the viscous time scale is of the order of 
microseconds. However, the sonic time scale is proportional to the size of the 
particle, whereas the viscous time scale is proportional to the square of the size 
of the particle, and already for a particle a few nanometres in size, the sonic 
time scale and the viscous time scale have the same order of magnitude - namely, 
several picoseconds. Clearly, the viscous time scale can further be reduced by 
employing a more viscous fluid, and conditions under which backtracking occurs 
can be reached experimentally. 

To complete the analysis, we study the velocity autocorrelation function of 
the solvent itself P^. For small time steps, it is known to exhibit anticorrela- 
tions, and it may be tempting to assume that anticorrelations in the velocity 
of the solute are a direct consequence of anticorrelations in the velocities of the 
solvent particles. However, we find that anticorrelations within the solvent and 
anticorrelations of the colloidal particle do not necessarily arise together, as can 
be seen by comparing the data for t//ts — 1.5 in [2] and |6] (left). 
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4 Conclusions 



We have confirmed that interference of sound and viscous effects leads to a non- 
monotonic decay of the velocity autocorrelation function of a single colloidal 
particle suspended in a viscous solvent. Furthermore, we have verified that 
under certain conditions backtracking can occur, i.e. the velocity autocorre- 
lation function becomes negative. This behavior is reminiscent of viscoclastic 
behavior - however, it has a very different, dynamical origin which lies in the 
interference of the formation of sound waves and hydrodynamic vorticity. Fi- 
nally, we have confirmed that despite the strong influence of backtracking on 
the velocity autocorrelation function, its long-time integral, the diffusion coef- 
ficient, is unaffected. That is to say that at long times the effects of sound 
are completely integrated out. An important consequence of this result is the 
confirmation that sonic effects at the single-particle level can be neglected at 
long observation times. 
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